Partial Kondo screening in frustrated Kondo lattice systems 
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We investigate the effect of geometrical frustration on the competition between the Kondo cou- 
pling and the Ruderman-Kittel-Kasuya-Yosida interaction in Kondo lattice systems. By variational 
Monte Carlo simulations, we reveal an emergent quantum phase with partial ordering in which the 
frustration is relieved by forming a magnetic order on a sublattice and leaving the rest in the Kondo 
screening with spin-singlet formation. The role of quantum fluctuations and spin-charge interplay 
is elucidated. 
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Kondo lattice systems provide a fertile ground for 
studying fascinating phenomena in strongly-correlated 
electron systems The key concept is competition 

between the Kondo coupling and the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction. The former is a 
local antiferromagnetic (AFM) interaction between con- 
duction electrons and localized moments, which promotes 
spin-singlet formation and results in Fermi liquid states, 
as the so-called Kondo effect . On the other hand, the 
latter RKKY interaction is an effective magnetic coupling 
between localized moments mediated by conduction elec- 
trons, which tends to stabilize a magnetic ordering Q. 
The competition leads to a quantum critical point (QCP) 
between a Fermi liquid state and a magnetically ordered 
state QCP has attracted much attention as a source 
of fascinating phenomena, such as non-Fcrmi-liquid be- 
havior and superconductivity. 

In the present study, we explore yet another phe- 
nomenon emergent from the competition between the 
Kondo coupling and the RKKY interaction. Our inter- 
est is in the possibility to have an intermediate quan- 
tum phase induced by geometrical frustration, which pre- 
empts QCP, with a coexistence of screened local moments 
due to the Kondo singlet formation and magnetically 
ordered moments stabilized by the RKKY interaction. 
This is a partially ordered state, which we call the par- 
tial Kondo screening (PKS) state in this paper. 

The partial ordering is sometimes seen in localized spin 
systems with geometrical frustration Q . Our target PKS 
is, however, qualitatively different from the partial order- 
ing in localized spin systems in the following points: The 
system includes itinerant electrons, and the disordered 
sites are not simply paramagnetic but participate in the 
Kondo singlet formation with vanishing their moments. 
These differences will bring about many distinctive as- 
pects not only in the stabilization mechanism of the par- 
tial order but also in the resulting physical properties. 

Several candidates for PKS are experimentally found 
in /-electron compounds, e.g., a distorted kagome mate- 
rial CePdAl and a triangular material UNi4B 
These PKS states were theoretically studied by the mean- 
field approximation of a pseudomoment model, which de- 
scribes the magnetic and singlet states by discrete classi- 
cal variables |10l-[l2|. In the previous studies, the effects 



of quantum fluctuations and the interplay between con- 
duction electrons and localized moments are not fully 
taken into account, despite the fact that they are obvi- 
ously crucial in the spin-charge coupled systems. 

In this Letter, we explore PKS as a quantum phase for 
the Kondo lattice model and the Kondo necklace model 
on frustrated lattices. The Kondo lattice model (KLM) is 
one of the fundamental models for rare-earth compounds, 
whose Hamiltonian reads 



{ij)tT i (ij) 

a) 

where the first term describes the hopping of conduction 
electrons and the second term represents the Kondo cou- 
pling between the conduction electron spin Ti and the 
localized spin Si. For simplicity, we consider S — 1/2 
spins for the localized spins. In Eq. ^ wc extend the 
model by adding the last term, the AFM Ising interac- 
tion between localized spins {S^ is the z component of 
Si), in order to mimic the magnetic anisotropy often seen 
in real materials [l^ . The Kondo necklace model (KNM) 
is a simplified variant of KLM at half filling Q : 
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where the charge degree of freedom of conduction elec- 
trons is suppressed with assuming that there is one elec- 
tron localized at every site. Instability toward PKS in 
KNM was studied previously by the authors for limited 
lattice geometries and system sizes [13] ■ We consider the 
two models on one of the simplest frustrated lattices, the 
two-dimensional triangular lattice, and take the sums {ij) 
over the nearest-neighbor sites. 

Wc study the ground state of the models ([T]) and ([2]) 
by the variational Monte Carlo (VMC) method. The 
method has several advantages compared to others; e.g., 
it takes account of quantum fluctuations neglected in the 
mean-fleld approximation, and it can avoid the minus 
sign problem even in frustrated systems. We here con- 
sider the variational wave function in the form 



(3) 



which describes magnetic states and nonmagnetic sin- 
glet states on an equal footing as a natural extension 
of the Yosida-type wave function [l5|. Here Impair) is 
a generalized BCS wave function defined by |</>pair) = 
iJ2e!L=i femal^al^^)^''\0), where fem are the variational 
parameters and |0) is a vacuum: The fermion operators 
aicr represent both conduction Cia and localized electrons 
(or localized spins), resulting in pair creations for any 
combination of them. We focus on the half-filling case by 
setting Np, — N, where N is the number of lattice sites. 
The electron number is fixed to be one at each localized 
spin site. C^^^ and C^^^ arc the quantum-number pro- 
jection operators for the total spin singlet and the total 
momentum zero, respectively. Pq is the Gutzwiller factor 
for optimizing the weight of configurations with double 
occupancies in KLM; Pq = exp{— J2i gin-ifn-ii,) where gi 
are the variational parameters and riicr = c[^Cia- We fol- 
low Ref. [l^ to optimize a large number of variational 
parameters by using the stochastic reconfiguration [l7j 
and to enforce the quantum-number projections. 

In the present study, we explore the solutions with 
three-sublattice ordering by imposing £^=° only for the 
same sublattices and by considering only the sublattice 
dependence of gi in Pq. C^^^ is used only for /, = 

0. Typically, the optimization is achieved by 300-1000 
stochastic reconfiguration steps with 1600-6000 Monte 
Carlo samplings. We confirm that the optimized wave 
function gives a precise ground-state energy compared to 
the results by the exact diagonalization; e.g., the relative 
error is less than 0.03 for KNM with N = 12. We apply 
the method to the clusters with system size N — 12, 18, 
and 24 with imposing the boundary conditions compati- 
ble with the three-sublattice order. 

First we discuss the results for KNM given by Eq. ^ . 
Figure [TJa) summarizes the ground-state phase diagram 
determined by VMC. There are three distinct regions, 

1. e., a magnetically-ordered (MO) state in the small J jW 
region, a Kondo spin liquid (KSL) state in the large J jW 
region, and a PKS state in between. The MO state for 
small J jW has a three-sublattice ordering as schemati- 
cally shown in the figure. This peculiar order is governed 
by the first and third terms in the Hamiltonian Q as 
discussed later. In the opposite large- J/P^ region, the 
second term becomes dominant, and the Kondo singlet 
is formed at every site to realize KSL. In the intermediate 
competing regime, we obtain the PKS phase in which one 
sublattice is dominated by the local Kondo singlet forma- 
tion and the remaining two retain magnetic ordering, as 
schematically depicted in the figure. 

Then we show how we identify three phases in the fol- 
lowing. Figure [2] plots J jW dependences of the onsitc 
correlation (x^ - Si) p_ and the intersite correlation between 
the localized spins on the same sublattice {Si ■ Sj)fj_ (/i 
denotes A, B, or C). In the intermediate J/W region, 
the onsite correlation (r^ • S^)^ on one sublattice (C in 
this case) becomes considerably larger in magnitude than 
the other two. as shown in Figs. [Ua)-(c). The difference 
A = \{n ■ Si)c - {{Ti ■ S^)a + in ■ S^)b)/2\ is plotted 
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FIG. 1: (color online). Phase diagrams for (a) the Kondo 
necklace model and (b) the Kondo lattice model at half filling 
on the triangular lattice determined by VMC. The lines are 
guides for the eye connecting the data for A*' — 18. There 
are three phases: MO, PKS, and KSL. Schematic pictures 
for the spin state in each phase are shown in (a), where the 
gray (black) arrows represent Ti (Si) and the circles denote 
the Kondo singlets. The dot-dashed and dashed lines in (a) 
show the phase boundaries for MO-PKS and PKS-KSL, re- 
spectively, determined by the mean-field approximation. 



for different system sizes A'^ in each inset. In the same 
region, {Si ■ Sj)c is suppressed compared to those for A 
and B sublattices, and moreover, it decays quicker with 
increasing distance than the other two [Figs. [Hd)-(f)]. 
These are clear indications of PKS; the local Kondo sin- 
glet is enhanced on the C sublattice compared to the rest 
magnetically active sublattices [l8| . 

Such specific spin configuration in the PKS state is 
also evidenced by the intersite correlations among differ- 
ent sublattices plotted in Figs. [3l^a) and[3l^b). For both 
Ti and Si, the intersite correlations measured from the 
C sublattice arc suppressed, while those between A and 
B sublattices arc robustly AFM. Therefore, the spin con- 
figuration of the PKS state is basically composed of en- 
hanced Kondo singlets on one sublattice and the AFM 
network on the remaining unfrustrated honeycomb lat- 
tice, as schematically shown in Fig. [IJa). 

In the smaller J/W region. Figs. [3{c) and Eljd) show 
that {Ti-Tj) is almost independent of the sublattice, while 
(S, • Sj) ~ ±1/4 between A and B but {Si ■ Sj) ~ be- 
tween C and other sublattices. These behaviors indicate 
the three-sublattice order for the MO phase in Fig.[TJa): 
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FIG. 2: (color online), (a)- (c) Onsite correlation for the three- 
sublattice ^ = A, B, C. The insets show their difference A 
for various system sizes A'^. (d)-(f) Intersite correlation be- 
tween localized spins on the same sublattice. Closed and open 
symbols show the results for the distance -\/3 and 3, respec- 
tively, as shown in the inset in (d). Data in the main panels 
are for KNM with iV = 18 and at (a),(d) = 0.0, (b),(d) 
l^jW = 0.4, and (c),(f) I^/W = 0.8. The lines are guides for 
the eye. 



Ti forms almost 120° order to optimize the W term, while 
Si exhibits an almost collincar AFM order in A and B 
sublattices with leaving perpendicular spins on the C sub- 
lattiee to optimize the Iz term. 

On the other hand, for large J/W, all (r^ • 5^^}^ become 
identical within the statistical errorbar and take a value 
close to —3/4, as shown in Figs. [5Ja)-(c). At the same 
time, all {Si • Sj)^ merge and become small with showing 
a rapid decay with the distance, as plotted in Figs. [2jd)- 
(f). These are a sign of the KSL phase where the spins 
form singlets locally as (xi ■ 5^)^ ~ —3/4. 

The phase boundaries are determined from the behav- 
iors of these spin correlations. In particular, the MO- 
PKS boundary is determined by a sudden change of in- 
tersite spin correlations in Figs.[lljd)-(f) as well as by the 
onset of A in Figs. [5l^a)-(c). Two states have the same 
symmetry in terms of spins apparently [l8| , but we con- 
clude that there is a phase transition between them in 
the light of the mean-field results discussed below. 

Let us make a remark on the limit of ^ J (3> W). 
In this limit, Si are decoupled from ri and expected to 
be disordered with macroscopic degeneracy [ly. Since 
J tends to lift the degeneracy through flipping Si, the 




FIG. 3: (color online). Intersite correlation at Iz/W = 0.4 for 
(a),(b) J/W = 0.6 and (c),(d) J/W = 0.2 in the KNM with 
N = 18. The correlations are measured along the horizontal 
arrows in the inset in (c) with fixing the site i at the origins 
of the arrows and shifting j. 



system is mapped onto a transverse-field Ising model on 
the triangular lattice, for which a three-sublatticc partial 
order is suggested to appear [23|. Such consideration 
leads us to expect the three-sublattice PKS state in the 
limit of Iz ^ J; this appears to be consistent with the 
result in Fig. [TJa). 

The VMC phase diagram is compared with that by the 
mean-field approximation (MFA) in Fig. [TJa). Here we 
perform MFA which accommodates three-sublattice or- 
dering by extending the method in Ref. MFA also 
predicts three phases, MO, PKS, and KSL. However, 
there are several differences compared with the VMC re- 
sults. The most distinct one is that VMC predicts the 
PKS state in a finite range of J/W down to Iz = 0, 
whereas it disappears at = in MFA. Another im- 
portant difference is that VMC phase boundaries shift 
to smaller J/W compared to the MFA results, and the 
width of MO state becomes relatively narrower. These 
are consequences of intersite quantum fluctuations ne- 
glected in MFA. 

Now we turn to KLM in Eq. ([1]). which includes itin- 
erancy of electrons explicitly. Figure [Ijb) summarizes 
the phase diagram at half filling. We find that KLM 
exhibits a similar sequence of three phases including 
PKS. The corresponding onsite correlations are plotted in 
Figs.|li;a)-(c). As seen in KNM [Figs. [i;a)-(c)], (r, • S^)^, 
becomes larger on one sublattice than the other two in 
the PKS region. Although A [insets in Fig.gjb) andgl^e)] 
suffers from finite-size effects (the reason is discussed be- 
low), we determine phase boundaries from qualitatively 
similar changes of spin correlations to those in Figs.[2{d)- 
(f) and in Fig. [3] (not shown). At Iz = 0, PKS behavior 
is not clearly observed in KLM [Fig. |4l[a)]. Furthermore, 
the MO region is wider than in KNM; we return to this 
point below. Consequently, PKS region becomes rela- 
tively narrower compared to the KNM case, yet remains 



4 



0.0 

-0.1 
3.-0.2 
of -0.3 
V -0.4 

-0.5 

-0.6 
0.0 







- A A > 




□ B 




- oC 




-(a) 





1.3 
1.2 

V 1.0 







A A \ 




□ B ^ 




oC ^ 














FIG. 4: (color online), (a)-(c) Onsite correlation and (d)-(f) 
local electron density in KLM with = 18. The insets show 
their differences with varying TV. The data are for (a),(d) 
h = 0.0, (b),(d) h/t = 0.4, and (c),(f) h/t = 0.8. 



robust between the MO and KSL phases even when con- 
duction electrons are considered explicitly. 

The relatively wider MO phase is presumably at- 
tributed to the complicated role of J in KLM: J enhances 
the RKKY interaction which tends to stabilize magnetic 
ordering, in addition to the enhancement of the spin- 
singlet formation. In addition, the notable system-size 
dependence in the small J /t region [the insets in Fig. |4] 
might be due to the long-ranged and oscillating nature 
of the RKKY interaction Q, which is difficult to cap- 
ture within the small finite-size clusters. Further study 
in larger system sizes is desired to clarify the nature of 
MO in KLM. 



Another interesting observation related with the 
charge degree of freedom in KLM is that PKS accompa- 
nies charge disproportionation. The local charge (ni)^ = 
'^^{nia) II disproportionates among the sublattices as 
shown in Figs. |4l^d)-(f). An instability toward charge 
density wave was recently discussed for the unfrustrated 
KLM model around quarter filling |2l|. The relation 
is not clear between the instability and our PKS with 
charge disproportionation. It is interesting to study the 
possibility of PKS for general filling, in particular, at 
commensurate filling. 

In summary, we have investigated the effect of geomet- 
rical frustration on the competition between the RKKY 
interaction and the Kondo coupling by the variational 
Monte Carlo simulation for the Kondo lattice model and 
the Kondo necklace model. The comparative study be- 
tween two models reveals the following features, (i) 
Both models exhibit PKS phase in between MO and 
KSL phases, (ii) The PKS state is further stabilized 
by quantum fluctuations and the spin anisotropy. (iii) 
Charge degree of freedom manifests in the stability of 
the MO phase and in the charge disproportionation as- 
sociated with PKS. All the results illuminate the appear- 
ance of the PKS phase in the Kondo lattice systems even 
when taking into account quantum fluctuations and spin- 
charge interplay, both of which have not been studied in 
the previous studies. We believe that the results pave the 
way for further understanding of PKS observed in com- 
plicated materials and of expected spin-charge entangled 
phenomena inherent to PKS. Our results will also cast a 
new light on the recent efforts to explore a new paradigm 
of QCP physics in the Kondo problem [2^ [23j. 
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